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We analyze effects due to transverse degrees of freedom in QCD calculations of the fundamental 

C^ ' hard exclusive amplitude of 7*7 -^ tt" transition. A detailed discussion is given of the relation 

^ . between the modified factorization approach (MFA) of Sterman et al. and standard factorization 

■ ' (SPA). Working in Feynman gauge, we construct basic building blocks of MFA from the one- 

_^ . loop coefficient function of the SFA, demonstrating that Sudakov effects are distinctly different 

^Nl ' from higher-twist corrections. We show also that the handbag-type diagram, contrary to naive 

c 2 ^ , expectations, does not contain an infinite chain of [M"^ /Q"^)" corrections: they come only from 

[^ ' diagrams with transverse gluons emitted from the hard propagator. A simpler picture emerges 

within the QCD sum rule approach: the sum over soft qG . . . Gq Fock components is dual to qq 

r^^ ' states generated by the local axial current. We combine the results based on QCD sum rules with 

Mh, pQCD radiative corrections and observe that the gap between our curves for the asymptotic and 

O ,' CZ distribution amplitudes is sufficiently large for an experimental discrimination between them. 
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I. INTRODUCTION 

The form factor F^.^*^o(q^, g^) relating two (in general, virtual) photons with the lightest hadron, the pion, plays 
a crucial role in the studies of exclusive processes in quantum chromodynamics. With only one hadron involved, 
it has the simplest structure analogous to that of the form factors of deep inelastic scattering. At large photon 
virtualities, comparing the pQCD predictions |p]-pl with experimental data, one can get important information 
about the shape of the pion distribution amplitude if nix). Due to its relation to axial anomaly 0], the 7*7* 7r° 
form factor was an object of intensive studies since the 60's P-jl^. Experimentally, F^-^,n°{lijl2) ^'^^ small 
virtuality of one of the photons, qf « 0, was measured only recently at e+e^ colliders by CELLO ||lj] and CLEO 
|L5[ collaborations (in the latter case, only a preliminary announcement of the results was made). The possibility 
to measure Fj-y-n" (9i ~ 0, g^) at fixed-target machines like CEBAF of Jefferson Lab was also discussed [M. These 
measurements inspired the studies of the momentum dependence of this form factor within various models of the 
nonperturbative quark dynamics p7|-|27[|. 

For a detailed comparison of pQCD predictions with experimental data, one should have reliable estimates of 
possible corrections to the lowest-order handbag contribution, in particular, those due to the gluon radiation and 
higher twist effects. Within the standard pQCD factorization approach, the one-loop radiative corrections to the 
coefficient function were calculated in refs. iH-n]. The authors of refs. [ pi|j28| incorporated the modified factorization 
approach of Sterman and collaborators [ p9| , p0[ in which the factorization formula invloves an extra integration over 
the impact parameter b± and Sudakov double logarithms of (a^ In (^j.))" type are summed to all orders. In 
refs. pl| , |2q | it was claimed that such an analysis takes into account some transverse-momentum effects neglected 
within the standard factorization approach p]j3l|-p3^. Incorporating the transverse- momentum-dependent wave 
function ^(x, fcj^), Jakob et al. ^l| also proposed a model for the effects due to the intrinsic (primordial) transverse 
momentum. 

Another attempt to take into account the transverse momentum effects was made by Cao et al. |P6| where the 
light-cone formalism expression M for the 7*7 -^ tt^ was used. Adopting an exponential ansatz for the transverse 
momentum dependence of the wave function, the authors observed large "higher- twist" corrections, with the 
conclusion that it is difficult in such a situation to make a clear distinction between different shapes of the pion 
distribution amplitude. 

In this paper, we will discuss various types of transverse momentum effects for the 7*777*^ form factor. First, 
we briefly outline the derivation of the leading-twist pQCD formula for this process using a covariant OPE-like 
factorization approach p4,pl^,B5[. In this framework, we identify the basic types of the higher twist corrections 
neglected in the leading-twist approximation. We show, in particular, that for massless quarks in a scalar theory 
no intrinsic transverse momentum effects are neglected in the handbag diagram: due to the simple singularity 
structure of the massless quark propagator, such effects can be taken into account exactly and lead to negligible 
pion mass corrections rn^/Q^ only. In QCD, the handbag diagram contains a twist-4 term interpretable as a O(fc^) 
correction, but no terms corresponding to higher powers of k'j^. Hence, the infinite tower of (M^/Q^)" corrections 
is generated by operators corresponding to higher qG . . .Gq Fock components. In Section II, we also discuss the 
structure of the results for the one-loop radiative corrections [Q-^ calculated within the standard factorization 
approach |l§,||,||,|3. 

In section III, we give a detailed one-loop derivation of the basic formulae of the modified factorization approach 
(MFA). We write the relevant one-loop integrals in Sudakov variables used in [ p9pC| ], introduce the impact pa- 
rameter b±_ as the Fourier conjugate variable to the transverse momentum fcj^ and reproduce (at one loop) the 
structure of the modified factorization |29]. In contrast (and complementary) to the original analysis, we use Feyn- 
man gauge which allows us to make a direct graph by graph comparison with the results ^^ obtained within the 
standard factorization approach (SFA). Since the modified factorization formulas appear as an intermediate step 
in our calculations which eventually produce the results of the SFA, the two types of factorization give identical 
results at any finite order of perturbation theory. The difference between the two approaches is only in different 
organization of all-order summation of higher- loop terms. Namely, in the MFA, the Sudakov- type double log- 
arithms (tts In {Qb±))"' are treated as logarithmic enhancements and are summed over all orders to produce a 
factor suppressing the contributions from the large-6 region. In the standard approach, the (a^ In (Q&))" terms 



are integrated over b± and included order by order. We show that for the 7*77r° form factor the use of the SFA 
procedure is well justified since the results of the 6 j_ -integration produce rather mild corrections ( ^ 20% at one 
loop). Another lesson from our detailed one-loop study of the MFA is that though the factorization formula of 
the MFA explicitly involves an integral over the impact parameter bj_ (or transverse momentum k±), the results of 
such an integration do not produce power suppressed contributions. Thus, despite the claims made, e.g., in refs. 



| p8 39 21I| higher-twist corrections are not included in the MFA . 

In section IV, we discuss two recent attempts pl| , p6[ t o model the intrinsic momentum corrections for the 
F^.^^o((5^) form factor. The approach of Jakob et al. [^l[ is based on the extrapolation of the modified factor- 
ization formula into the nonperturbative region. At large impact parameters 6, the Sudakov suppression factor is 
supplemented by the nonperturbative wave function ^{x,h) reflecting the effects due to the primordial transverse 
momentum distribution. However, since terms which were inessential for the derivation of the Sudakov factor at 
large Q^ may be quite important for small Q^ , it is not clear for which Q^-region such an extrapolation is sufficiently 
accurate. We observe, in particular, that instead of producing the Q^ = value dictated by the axial anomaly 
MG], the extrapolation formula gives a logarithmically divergent result suggesting that the extrapolation should 



not go down to very low Q^ . Cao et al. |g6| use the expression for the qq Fock state contribution to F^<.^j^o{Q'^) 
derived in the light-cone formalism by Brodsky and Lepage ||I|. This expression involves no approximations and 
has correct limits both for small and large Q^. In particular, we demonstrate that, in full accordance with our 
general analysis, it contains no higher-twist contributions. Still, one should take into account that the qq term, by 
definition, does not include the contribution due to higher qG . . .Gq Fock components of the pion light cone wave 
function. As shown in ref. g , the latter coincides in the real photon limit Q^ = with that of the qq Fock com- 
ponent and doubles the total result at this point. Clearly, the inclusion (or at least modelling) of this contribution 
is necessary for a consistent description of subasymptotic effects. Comparing the approaches of refs. |21^,Eq], we 
emphasize that they incorporate two completely different light-cone schemes. The light-cone formalism of Brodsky 
and Lepage used in ref. pq ] is equivalent to incorporating the infinite momentum frame. On the other hand, the 
approach of ref. [|lj (and that of the underlying papers ^,Q ) is based on the Sudakov decomposition. The basic 
difference between the two light-cone approaches is that the momentum of the virtual photon in the 7*7 — > tt^ 
process is dominated by the transverse component in the BL light-cone scheme while it is purely longitudinal in 
the Sudakov approach. 

In Section V, we use QCD sum rule ideas to get a model for the F^.^^o((3^) form factor which reproduces 
both the Q^ ~ constraint imposed by the axial anomaly and the lowest-order pQCD results for high Q^. We 
show also that the results obtained on the basis of QCD sum rules and quark-hadron duality can be interpreted 
in terms of the effective valence wave function which absorbs information about soft dynamics of higher Fock 
components of the standard light-cone approach. Combining these results with pQCD radiative corrections, we 
obtain an expression depending on the choice of the low-energy distribution amplitude. The difference between our 
results for the asymptotic and CZ distribution amplitudes is sufficiently large for an unambiguous experimental 
discrimination between these two possibilities. 

II. FACTORIZATION 
A. Structure of factorization 

We define the form factor F^<-y*-^o (^qf, q'^^ of the 7*7* -^ 7r° transition through the matrix element 

An J {7r,P \T{J^iX)MO)}\0)e-'^'''d^X = ie^V2e^,a.pq1q^F^,^,,o (qlql) (2.1) 

where J^ is the electromagnetic current of the light quarks 

J^L = e„U7^u -I- edd-ff^d (2.2) 

and \n, p) is a one- pion state with the 4- momentum p. Note, that our definition (aimed at getting a simple 
coefficient for the spectral density for the triangle anomaly diagram, see Section V) differs from that in refs. 



pl|j2q| by factor ■\/2/47r. Experimentally, the most favorable situation is when one of the photons is real or 



almost real: qf ^ 0. In this case, we will denote the form factor by F^,^^o{Q^), where Q^ = — gf is the virtuality 
of the other photon. It should be sufficiently large for pQCD to be applicable. In general, a power-like behavior of 
F^*^Tro{Q'^) in the large-Q^ limit can be generated by three basic regimes (see Fig.l). 




q^^^[yy0^ p q^ 
b) 




c) 



FIG. 1. Structure of factorization for the F^,^^o{Q'^) form factor at large Q^. 

The dominant contribution is provided by the first regime (Fig. la) which corresponds to large virtuality flow 
through a subgraph V containing both photon vertices. The power counting estimate for the large-Q^ behavior of 
such a configuration with arbitrary number of external lines of V is given by (see refs. [p2|,p4[) 



FiQ') < Q 



■E' 



(2.3) 



where t^'s are twists (dimension minus spin) of the quark and gluon external lines of V, with t — 1 for the quarks 
and t = for the gluons in a covariant gauge. Hence, for the leading term, one should take the minimal number 
of quark lines (two in our case) while the number of the gluonic ^-fields is arbitrary. Generically, the leading 
contribution of this type can be written as 



F^,^,^oiq,,q2) = J Ci^,r^,ql,q2■,^i^){p\OiC,v)m^.-d^Cd^rJ, 



(2.4) 



where the parameter /i^ is the factorization scale, (7(^,77,(71,(72) corresponds to the short-distance amplitude with 
two external quark lines and 0(^,77) is a composite operator 0{^,rj) ^ q{(^)'^5^^E{^,rj; A)q{ri). The path-ordered 
exponential 



E{^,j];A) = Peii-piig A^{z)dz^' 



of the gluonic field A results from summation over external gluon lines of V. For the quark propagator, e.g., one 
has 

S'{^ -V)+Js^i^- z)rgA^{z)S-{z -7j)d^z + ...^ E{^, 77; A) S^i^ - 77) [l + 0{G)] (2.5) 

where 0{G) depends on the gluonic fields through the gluon field strength tensor G^jy and its covariant derivatives. 
Since G^i, is asymmetric with respect to the interchange of the indices /^, i/, it should be treated as a twist-1 field. 
Basically, the contribution ( ^.4| ) is analogous to the quark-antiquark term of the standard operator product 
expansion for J"(0)J^(z). In this form, the operator 0(^,77) still contains non-leading twist terms. To get the 
lowest-twist part, we should expand 0{£,, rj) into the Taylor series 



00 _. 
mi^l^m, W, A)q{ii) = ^ -A'-i A-^ . . . A''-q{0l5luD,,D,, . . . D,^q{0 ; A = 77 - C 



(2.6) 



n=0 



and pick out only the symmetric-traceless part qjzijjjD^-^D,^^ ■ ■ ■ -Diy„}(? of each local operator from this expansion. 
The traces correspond to operators with contracted covariant derivatives D^ D^ which, for dimensional reasons. 



are accompanied by powers of the interval (^ — 77)^. Likewise, the (^ — 77)^ factors produce extra powers of z^ after 
integration over ^ and r\. Finally, each power of z^ results in an extra power of l/Q"^ , i.e., each pair of contracted 
covariant derivatives D" ...D^ in a higher-twist operator produces 1/Q^ suppression at large Q^. Hence, the 
twist-2 part of 0(1^, 77) corresponds to the lowest term of the expansion over (1^ — rj)^ 

0(e,r;) - 0(e, 77)|(5_,).^o + 0((e - V)^)- (2.7) 

The light-cone matrix element can be parametrized in terms of the pion distribution amplitude (DA) ifj^^x): 



(oia(e,r?)k",p) 



(e-,,)2=o 



ip^ I e 



-ix{^p)—ix{rjp) 



ip^{x)dx , 



(2., 



which gives the probability amplitude that the fast-moving pion is a qq pair with its longitudinal momentum p 
shared among the quarks in fractions x and x = {I — x) (throughout the paper, we use the "bar" convention for 
the momentum fractions: x=l — x,y=l — y, etc.). Substituting this representation into the generic expression 
( ^.4| ), we obtain the hard scattering formula 



i^. 0(51,(72) 



47r 



Tiqi, q2; xp, xp) (p^{x) dx 



(2.9) 



where the factor 47r/3 is due to our normalization of the form factor and T{qi,q2;k,k) is the amplitude for the 
subprocess l{qi)l*{q2) -^ 'l(k)q{k). Calculating this lowest-twist amplitude in the momentum representation, we 
should realize that the neglect of the higher-twist operators having extra D^ is equivalent to taking fc^ = 0, fc^ = 
for the external quark momenta. In general, this limit is singular for diagrams with loops, and one should regulate 
the resulting mass singularities Infc^ in some way, e.g., by dimensional regularization or by taking massive quarks 



and k 



In the latter case, only the logarithmic r7T,g-dependence should be kept in the final result: keeping the 



power terms rnz/Q exceeds, for light quarks, the accuracy of the method. The subsequent procedure is to split 
the logarithms \n{Q^/m^) into the long-distance and short-distance parts ln{Q^/m^) — ln((5^//i^) + ln(/i^/7n^) and 
absorb the long-distance ones ln(/i^/TO^) into the pion distribution amplitude: ipT^ix) —> ip.n{x;iJL). 



^1 



li 



(1-x)p 



xp 



FIG. 2. Lowest-order diagram. 

Thus, the lowest- twist contribution corresponds to the parton picture in which only the longitudinal (proportional 
to p) components of the external quark momenta appear. In the lowest order (see Fig. 2), the amplitude for 
transition of two photons into the quark-antiquark pair with coUinear lightlike momenta xp, xp is givenj] by the 
quark propagator: 



To(x,Q2 



1 



1 



-((71 - xpY xQ^ 



(2.10) 



^In fact, there are two diagrams obtained from one another by the interchange of photon vertices. However, due to the 
symmetry of the distribution amplitude ipTr{x) = <^^(1 — a;), their contributions can be united. 



and the pQCD result ^ for the large-Q^ behavior of the form factor is 

F,',AQ ) = ^J^-^da: = -^lo. (2.11) 

Necessary nonperturbative information is accumulated in the same integral 



Io = ^ [ ^^ dx=^ f To(x, Q2) ^^(a;) dx (2.12) 



that appears in the one-gluon-exchange diagram for the pion electromagnetic form factor [p3,BllE2J. The value of 

nj4T 



I depends on the shape of the pion distribution amplitude ipTr{x). In particular, using the asymptotic form 32 1 

if^^'ix) = 6Ux{l - x) (2.13) 

gives Iq'^ = 3. If one takes the Chernyak-Zhitnitsky ansatz B^] 

ip'^^ix) = 30Uxil-x){l-2x)\ (2.14) 

the integral /q increases by a sizable factor of 5/3: Iq^ — 5 and one can hope that this difference can be used for 
an experimental discrimination between the two competing models for the pion DA. 

Since one of the photons has a small virtuality, one should, in principle, also take into account the regime (see 
Fig. 16) involving a long-distance propagation in the gi-channel, with large momentum flowing through a central 
subgraph W containing only the virtual photon vertex. In the lowest order, this subgraph corresponds to a hard- 
gluon exchange, just like in the asymptotically leading pQCD contribution to the pion electromagnetic form factor. 
The power counting for such a contribution into F^,yT^o[Q'^) is given by 

^(Q^) ^<9"*°i~*°^ , (2.15) 

where t^. , i = 1, 2 are the twists of composite operators Oj corresponding to qi- and p-channel, respectively. Taking 
into account that twist of a gauge-invariant color-singlet composite operator Oi cannot be less than 2, we conclude 
that this regime gives a nonleading 0{1/Q*) contribution. 

The third regime (Fig.lc) corresponds to Feynman mechanism, i.e. to a situation when the passive quark is soft. 
Using the wave function terminology, we can say that F^t^^oiCp) in this regime is given by an overlap of soft wave 
functions describing the initial and final state. This contribution also behaves like l/Q^ at large Q^ . 

B. Handbag diagram and transverse momentum 

For the OPE contribution, the simplest power corrections come either from the traces of the two-body operator 
C(a;, y) which appears in the handbag diagram or from a direct insertion of gluon lines with physical polarizations 
into the propagator connecting the photon vertices. Since D'^ Dn can be interpreted in the momentum representation 
as the (generalized) virtuality k"^ of the quark field, the higher-twist operators containing D'^D^, look like a natural 
candidate for description of the effects due to the transverse momentum of the quarks. However, there are some 
practically important amplitudes which, due to their simple singularity structure, are "protected" from the towers 
of (D^)"-type higher-twist corrections. The most well-known example is given by the classic "handbag" diagram 
for deep inelastic scattering. The lowest-order diagram for the 7*7 -^ 7r° form factor (Fig. 2) has similar properties. 
Consider its analog in a toy scalar model 

F{q2,p) ^^J e-*«^^(O|0(O)</>(z)b)^ . (2.16) 

The first term in the 2:^-expansion for the matrix element 

(O|0(O)0(z)|p) = U^p) + z^Uzp) + [z^fUzp) + . . . (2.17) 

6 



corresponds to the twist-2 distribution amplitude while subsequent terms correspond to operators containing an 
increasing number of 9^'s. It is straightforward to observe that, while the twist-2 term produces the 1/Q^ contri- 
bution, the twist-4 term is accompanied by an extra z^-factor which completely kills the 1/z^-singularity of the 
quark propagator, and d'^z integration gives S'^{q — xp), which is invisible for large Q^ . The same is evidently true 
for all the terms accompanied by higher powers of z^ . This means that the handbag diagram contains only one 
term with a powerlike behavior for large Q^: it cannot generate higher powers of 1/Q^ which one could interpret 
as the ((fc^)/Q^)" expansion. Since only the z^ = projection of the bilocal operator survives, we can parametrize 



(O|(/.(O)0(z)b) = / </7(x)e-"(^P)da; + ... , (2.18) 

Jo 

where the dots stand for terms producing the "invisible" contributions, and write the lowest-order term as 

n'?.-p)--rr^^^-=-rr^^^--r^5^!^rf- (2.19) 

Jo {<l2 - xpY Jo [Qi - x-pY Jq xQ-^ -f xxp-^ 

Hence, the handbag contribution in this case contains only the hadron-mass corrections (cf. Q]), but it gives no 
information about finite-size effects. In the momentum representation, the origin of this phenomenon can be traced 
to the fact that a straightforward expansion of the propagator is just in terms of traceless combinations: 

1 '^^ on ^^^ on 

'^(i^i < i9i) E T^s^^'' ■ ■ ■ '^''"i^^- • • ■ ^mJ + (^m > \q\) E im^'i'' ■ ■ ■ '?'"{^'^^ ■ ■ ■ ^'^"^■ 



il-kY '" '^"t^O^l'y^' ■ -'-•■— -'---^^Jfc2). 



(2.20) 



The handbag contribution corresponds to |A;| < \q\, and this part of Eq.(2.20) without any approximation produces 
an expression equivalent to treating the fc-momentum as purely longitudinal k = xp. 

It is worth noting here that though the hadron-mass corrections have a powerlike behavior {p'^/Q^Y^, they 
should not be classified as higher-twist corrections: they result from the kinematic hadron-mass dependence of 
the lowest-twist contribution. For deep inelastic scattering, the possibility to calculate the target-mass corrections 
within the lowest-twist contribution is known as the ^-scaling phenomenon ||45| , [46| . As emphasized by K.Ellis et 
al. [H, the ^-scaling phenomenon can be also understood in terms of the primordial transverse momentmn, if 
one takes into account that, for the lowest-twist term, the transverse momentum distribution is totally due to the 
non-zero hadron mass, i.e., it has a purely kinematic nature and for this reason can be calculated exactly. The 
quark propagator in QCD has a stronger singularity z/z"^. As a result, the handbag-type contribution in QCD 
contains a twist-4 operator with extra D^ |47[ , but no operators with higher powers of D^ . 



One may argue that there is another part in Eq.(2.20), when k is large (i.e., \k\ > l^j). In this case, the fc-line 



corresponds to high virtualities. If such a large momentum goes directly into the soft hadronic wave function, 
the Q^-behavior of such a contribution repeats the fc-^-dependence of the soft wave function, i.e., very rapidly 
(say, exponentially) decreases with Q^ (see Section IVC below for an explicit illustration). A more favorable 
possibility is when the large momentum by-passes the wave function. Such a configuration can give a leading- 
power contribution. In the latter case, the large virtuality flows through several lines forming a subgraph with the 
same (minimal possible) number of external quark lines as the lowest-order leading twist contribution. In the QCD 
factorization scheme, the relevant contribution produces a part of a higher-order coefficient function (see Fig. la). 

C. One-loop radiative correction to the coefficient function 

At one loop, the coefficient function for the 7*7 -^ 7r° form factor was calculated in refs. |4|-||: 



T(.,Q-M^)^^{l + C.g 



3 1 \ 1 /^2 I 2\ 1 1 2 x\nx 9 



(2.21) 



(1-x)p 




, r^ 



d) 




b) 



e) 



c) 




f) 



FIG. 3. One-loop diagrams. 

In full compliance with the factorization theorems |gl|,|| (see also P4| , p^y49| ) , the one- loop contribution contains 
no Sudakov double logarithms In'^ Q^ of the large momentum transfer Q. Physically, this result is due to the color 
neutrality of the pion. In the axial gauge, the Sudakov double logarithms appear in the box diagram 3c but they are 
cancelled by similar terms from the quark self-energy corrections 3d, e. In Feynman gauge, the double logarithms 
In Q^ simply do not appear in any one-loop diagram. It is easy to check that the term containing the logarithm 
ln((3^//j,^) has the form of convolution 



"TTTT Cprr^i - + In a; 



27r V2 



1 



,V{tx)d^ 



(2.22) 



of the lowest-order ("Born") term T'o(^, Q^) = 1/CQ^ ^^^d the kernel 



Vi^,x) 



as 

2^ 



Cf 



e 



9{^<x)(l + 



C 



+ ^9{^>x)[l + 



^-x 



(2.23) 



governing the evolution of the pion distribution amplitude. The "-I-" -operation is defined here, as usual 



by 



[F(e, x)]+ = Fit x) - <5(C ~x)J FiC x) dC . 





(2.24) 



Since the asymptotic distribution amplitude is the eigenfunction of the evolution kernel V{tx) corresponding 
to zero eigenvalue, 



V{tx)ip'''{x)dx = 0, 



(2.25) 



the coefficient | -I- Inx of the ln(Q^/yU^) term vanishes after the x- integration with (/?°^(a;). Hence, the size of the 
one- loop correction for the asymptotic DA is /.t-independent and determined by the remaining terms. The /-integral 



Q2 p1 

I=— T{x,Q'^)ipTrix)dx 

JTT JO 



(2.26) 



(cf. Eq.(2.12)) then can be written as 



/I 



3 1 



5 a^ 
2^^ 2^ 



(2.27) 



-7/4 generated by the i In 



9/2 (see Eq.(2.21)) partiaUy compensated 
contribution and +1/4 



The negative coefficient —5/2 here conies from the constant term 

by two logarithmic terms which give together +2, with 

by — a;lnx/[2(l — x)] term. With Cp = 4/3, the net factor is [1 — |as/7r]. Hence, for cCs/tt « 0.1, the one-loop 

correction is less than 20% and the as/ir expansion looks "reasonably convergent" . Taking the CZ form for 1^9(2;; /i), 

we get 



\^{x,^)^ip^^ {x) 



,_C,g(^MgW) + | 



(2.28) 



Again, the negative coefficient —49/72 comes from the —9/2 term compensated by an increased contribution from 

the logarithmic terms: ^ In^ x gives +263/72 and — xlnx/[2(l — x)] gives 1/6. For fi — Q, the one-loop modified 

factor is [1 — ^(cks/ti")], i-e., the total correction is smaller than that for the asymptotic DA. Since the result is 

49 
//-dependent in this case, by an appropriate choice of /i, namely, taking /i — e'^™Q ~ 1.5 Q we can formally get a 

vanishing 0{as) correction. Then the one- loop expression for the form factor would coincide with the lowest-order 

formula, but with the distribution amphtude (p^^{x;fj,) evolved to the scale /x « 1.5 Q. However, at this scale, 

tpTr{x\ n) does not necessarily have the CZ form. To treat the evolution in a consistent way, we set the boundary 

condition that 93^^ (x; jjl) has the canonical CZ form ^p'^^{x) = 30 fn xx (1 — 2a;)^ at some specific scale fj, — Qo (the 

original derivation p3[ assumes Qo — 0.5 GeX^). Taking into account that ip^^{x) is a combination of two lowest 

eigenf unctions of the evolution kernel, we can write the solution of the evolution equation in the leading logarithm 

approximation: 



^^^(rr; M) = <^(a:) + K^(x) - <^(x)} 



lnQ2/A2 



21 72 //3o 



In/i^/A^ 



(2.29) 



where 72 = 50/9 is the relevant anomalous dimension and /3o = 11 — iNf is the lowest coefficient of the QCD 



/3-function. In what follows, we take Nf 

5a^1 

3 TT I 



3 and /3o = 9. Choosing fi — Q, we get for the /-integral (cf. ]5l|]) 



^lv^(^,Qo)=v?^(a;) -3p- 



1- 



InQg/A 



lnQ7A2 



50/81^ 



5 1- 



49 as 
lOSV 



lnQ2/A2 



lnQVA2 



50/81 



(2.30) 



Note that the In x term generates a larger positive contribution for ip^^{x) because ip'^^ (x) is more concentrated 
in the end-point region x ^ than ip'^^{x). Furthermore, if the distribution amplitude is extremely concentrated in 
the end-point region x ~ 0, a positive contribution from the i In x term dominates the correction and generates 
a large positive net effect. In such a situation, the one- loop correction vanishes only ii jj. ~ aQ with a < 1. 
The broader the DA, the smaller should be the parameter a which reduces the one-loop expression to the lowest- 
order one. Since the effective normalization scale is smaller for a broader DA, perturbative QCD applicability 
is postponed to higher Q^. One may speculate that this phenomenon simply indicates that for a broad DA the 
quark virtuality xQ^ is a more natural choice for the effective factorization scale than the photon virtuality Q^ 
(i.e., a ~ (x)) and pQCD is applicable only if the average xQ^ rather than Q^ itself is large enough. One faces a 
similar situation studying the pQCD contribution to the pion form factor. The average virtuality {xyQ'^) of the 
exchanged gluon in that case is essentially smaller than Q^ and one may question both the self- consistency and 
reliability of the pQCD analysis at accessible energies |52|j53tl . In ref . ^^ , it was argued that due to the Sudakov 
effects in the impact parameter space, the pQCD treatment of the lowest-twist one-gluon-exchange term for the 
pion form factor is self-consistenM at smaller Q^ than suggested by the estimates of the magnitude of the average 



*Note, that self-consistency of the pQCD expansion (small q^ corrections) for the lowest-twist term does not necessarily 
mean that pQCD is reliable, since power corrections (M^/Q^)" can still be large (see discussion at the end of Section V). 
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gluon virtuality xyQ^. One may expect that similar effects m.anifest themselves also in the 7*77r*' form factor. 
Indeed, our numerical analysis of the one- loop correction shows that taking a = 1 (rather than a — (x)) provides 
a good choice for the factorization scale. It is accompanied by a small one-loop correction even for a broad DA of 
CZ type. 

It is worth noting here that, even without incorporating the impact parameter representation, one can observe 
some traces of the Sudakov effects in the structure of the one-loop coefficient function in the region of small fractions 
X. As explained earlier, the one-loop term is obtained by calculating the 7*7 —f qq amplitude for massive on-shell 
quarks with subsequent absorption of the mass logarithms in the form ln.{fj,'^ / m'^) into the distribution amplitude. 
When the virtuality xQ^ of the quark line connecting the photon vertices becomes small, the vertex correction for 
the virtual photon (Fig. 3a) is dominated (in Feynman gauge) by the ojf-shell Sudakov double logarithm which can 
be written as 



27r 



CfIu- 



Sln.<?^ 



xQ^ 



where xQ^ is the virtuality of the hard quark. Of course, since this virtuality is parametrically of the order of Q^ 



we get only a single logarithm with respect to Q^, namely, {as/2Tr) Cp ln((5^/m^) In a: (cf. (|2.2l| )), just as required 
for factorization. However, if we write the sum of two terms 






In X 



91 ^'1 ■ 
2 in — ^ in X 



which dominate the small- a; region as 



47r 



Cf 



\u' 



-In 



Q 



21 



we see that it converts into the standard on-shell Sudakov double logarithm 

47r ni'^ 



when xQ^ 



Of course, the region where xQ^ is parametrically of the order of the IR cut-off m"^ is outside 



the formal applicability region of the factorization approach, and there is no surprise that double logarithms of Q^ 
appear there. Note the well-known difference as/27r — > as/47r between the off- and on-shell forms of the double 
logarithms. In higher orders, Sudakov logarithms are expected to exponentiate producing the Sudakov form factor^ 
exp[— (as/47r) Cp ln^((5^/m^)], and the region of very small xQ"^ is relatively suppressed due to Sudakov effects. 
This also means that taking /x^ ~ xQ^ in Eq.(2.21) is not an optimal choice, since it is accompanied by a negative 



rather than vanishing correction. Indeed, the original motivation to take a lower scale fj, < Q was to compensate 



the positive contribution from the In x term. However, taking /f^ ~ xQ^ in Eq.(2.21) for a wide DA generates a 
negative (— In x) term which over-kills the original positive ^ In x term and converts its sign in the net result. A 
negative correction, in its turn, suggests that a larger factorization scale is a better choice. This indicates that, for 
a broad DA, the typical distances probed in the hard subprocess are larger than those corresponding to 1/Q^ but 
smaller than those corresponding to the inverse of the average quark virtuality xQ^ . 

As we will see in the next section, the modified factorization pG] is similar to the choice ff^ ^ xQ^ and for this 
reason it is accompanied by a negative correction. We will also explicitly show that the latter, in full accordance 
with the MFA analysis 129], can be explained by Sudakov effects in the impact parameter space. 



^For the pion EM form factor, exponentiation of a similar combination (CFaa/47r)[In^(a;j/Q^/77i'^) 
in ref . p& was verified by a two-loop calculation Im . 



'ln^(Q^/m^)] suggested 
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III. ONE-LOOP RADIATIVE CORRECTIONS AND TRANSVERSE MOMENTUM 
A. Vertex correction for virtual photon and Sudakov effects 

To establish the connection between standard and modified factorization approaches, we give below a rather 
detailed discussion of the structure of the one-loop coefficient function using the Sudakov decomposition for the 
loop momenta. We use the same definition of transverse momentum k^ as in ref. |^,|0[, introduce the impact 
parameter 5_l and then translate our results into the 5i^-space. To be able to make a diagram by diagram comparison 
with ref. [g|, we use Feynman gauge. This also allows us to give an independent one-loop derivation of the fej^-space 
Sudakov effects which complements the general approach |2^] based on the analysis in the axial gaugeQ. We 
find it also instructive to demonstrate how the 6^-space double logarithms appear in a situation in which double 
logarithms of Q^ are absent in any diagram. 

We start with the diagram 3a which is the most natural suspect in a search for Sudakov effects in Feynman 
gauge. According to general rules, calculating the coefficient function one should assume that external quarks carry 
purely longitudinal lightlike momenta xp and xp. Using p and qi (abbreviated in this section to q for convenience) 
as the basic Sudakov light-cone variables, we write the momentum k of the emitted gluon as 

k^{S,-x)p + riq + ki^ (3.1) 

and then take the 77-integral by residue. After that, the contribution of Fig. 3a (and any other one-loop diagram) 
can be schematically written as 

t!'\x,Q') = ^Cf dd AMx,Q';^,k^)—^ =^CFt,ix,Q'). (3.2) 



2tt Jo ^ J '' '^ '^' ^' 2tt 2tt 

The internal amplitude Ma{x, Q^; ^, /cj_) for the diagram 3a is given by 



The fc_L-integral diverges both in the fcj^ ^ cxi and fcj^ — > limits. The ultraviolet large-Zcj, divergences (they are 
actually irrelevant to our analysis) are removed by the -R-operation, while the 1ow-A;_l collinear divergences can be 
regulated by taking massive quarks. In that case, k\ -^ k^ + rn? and the small-fcj_ divergence (collinear singularity) 
is converted into the mass logarithm \n.(Q^ /m^) generating the evolution of the pion distribution amplitude. The 
Sudakov effects are also related to the l/k\ singularity. It is easy to check that the coefficient in front of \/k\ in 
the singular part 

Uf'i..Cf^W - __i^_^_Ol_g (1) «K > ,) ,3^4) 

has the form of the product of the Born term 1/^Q^ and the relevant part 



of the evolution kernel ( p.23| ). Note, that calculating the evolution logarithm hiQ^ /rn? from cPkj_/k^, one can 
take k± — ("neglect k±") in all other places, in particular, in the denominator factor S^Q^ + k\/S^. However, 
nothing prevents us from going beyond the leading logarithm approximation. Keeping the /c^-terms, we can take 



** In a recent paper p5|, Li gave a covariant gauge derivation of the modified factorization for inclusive processes and 
heavy-quark decays. However, in technical implementation, his approach is quite different from ours. 
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into account those contributions which do not have logarithmic behavior with respect to rri^ or Q^. We will see that 
among them, there are "Sudakov" terms with a specific double-logarithmic dependence on the impact parameter 
b±, the variable which is Fourier-conjugate to the transverse momentum k±. To separate the contributions related 
to the evolution kernel from those corresponding to Sudakov effects, we first make the decomposition 



1 



1 



1 



1 



(3.6) 



[^Q' + fci/e] ^Q^ \iQ' + fci/e xQy (e - x)Q2 + kl/^ 

and notice that the denominator factor ^Q^ + fc^/^ reduces to xQ^ when ^ — x and k± = 0. Hence, we can write 



2tt 



tr'{x,Q' 



d^ / d^k_i_ 



£,/x 



kl [CQ' + fci/C] xQ'' 



d£. 



llx 



iO" ^k\H\k\{{i- x)Q^ + k\l^ 



llx 



(C - .t)Q2 + llll 



- 8{^-x)b\ki_) \ dC / d^k±- 



To disentangle the product of the delta- functions in ^ and fcj^ variables, we rewrite Eg. (p. 71) as 



(3.7) 



d^ 



d^k^ 1 / {l/x)e{^>x) 

iQ' + kl/i \kl\{i- x)Q^ + kill, ^ 



^'^'-'^'A knc-.w.km -''''' 



d^k± 



ki 



{C~x)Q^ + kl/C 



rfC 



where the combination 

{l/x)ei(>x) 



(e - x)Q^ + kl/i) , (e - x)Q2 + ki/i 



{i/x) 0{^ > x) f iC/x)9{C>x) 



ic-x)Q^ + kyc 



(3.8) 



(3.9) 



is an analog of the "plus" operation for the case when the transverse momentum is present. Similarly, the expression 

1 „o,. . /■ d^k. 



kl [(C - 2:)Q2 + kl/C\ 



5\k^) 



kl 



{C-x)Q^ + kl/C 



(3.10) 



can be interpreted as a "plus" distribution with respect to fc^. Extracting the pure 1/fc^-singularity from the 
(• • •)+ term in Eq.(|3~ 



{^/x)e{i>x) 



1 ({£,lx)e{i>x) 



kl I [(e - x)Q^ + fci/e] / ^ Q^ki V {i-x) ; + Q^Ui-x) [i{(, - x)Q^ + kl] 



1 



(e/s) e{i > x) 



we can write (|3.8D in the impact parameter representation as 



1 
f--9{x,Q^) = l-j dij B{^;bQ) 



K(e, x)L{bm) + Ea{x, i- bQ) + 5{i - x)Saix, bQ) 



d^b±. 



(3.11) 



(3.12) 



The function B{^; bQ) gives the Born term in the 6-space 



B{^;bQ) = 



-ik\b_\ 



27tJ ^Q^ + kl/^ 



d'k^=^K„[bQJ^a, 



(3.13) 
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where b — \b±\ and Kq{z) is the modified Bessel function. By L{bm) we denote a reguiarized version of the integral 



resulting from the first term in Eq.(3.11): 



1 /" ,T, e 



''.k \ h \ 



In particular, if the integral is regulated by l/k\ — + l/(fc5_ + "m?), then L{bm) ~ Ko{bm). The function L{bm) 
contains the mass logarithm ln(m&) multiplied by the relevant part Va{^,x) of the evolution kernel. As discussed 
in the preceding section, the mass singularity In(TO) must be absorbed (in the form ln(m//i), where /i is the 
factorization scale) into the redefinition of the distribution amplitude: ifTr{x) — > ipTr{x;fi). The second term in 
Eq.(^.ll) is given by the function E{x,^; bQ) which also contains the evolution kernel Va{S,, x) 



E.i..,m - -rj e-- ( ,_--- ,,, j/^. - - 



e Oi^>x) 



Ko{bQJi^-x)^ 



(3.15) 



+ 



It is easy to notice that both the Born term B{£^;bQ) and the evolution-related terms L{bm) and Ea{x,^;bQ) 
exponentially decrease at large 5, since the function Kg {b . . .) behaves like exp(— 6 . . .) in this limit. On the other 
hand, the "Sudakov" term 



accompanied by S{£, — x) in Eq.(3.12) has a completely different behavior at large b. Indeed, changing the variable 
( in the above integral as 1 — C = yS, we rewrite Eq. dS.lq ) in the form 



Saix;Qb) = — / d^k^ -2 / _ , ,2 ,-,^, = s{xQb). (3.17) 



-Jd k^ ^ J y-^^^/-2Q2 



According to this representation, the function s{xQb) vanishes as 6 — > 0. In the opposite limit of large impact 
parameters, it has a double-logarithmic dependence on b. To see this, we integrate first over y and then over k± 
taking into account that the factor (e*''"^''^ — 1) provides, in the limit of large b, an effective IR cut-off at k± ^ 1/6. 
As a result, we obtain the large-6 behavior of s{xQb) [g9| 

s(xQb)^^Jd'k^^^^^^^^ln(^^^^ J^lnl^^^^-^ln^xQb), 1/Aqcz5 » 6 » l/Q. (3.18) 

1/6 

To be on safe side, we included the 1/A^qcd 3> b restriction to emphasize that these results are only valid in the 
region where one can trust pQCD expressions for quark and gluon propagators. 

Integrating s{xQb) with the Born term gives, for small x, a negative double logarithm — ^In x. As discussed 
above, such a correction is expected when one uses xQ^ as the factorization scale. Indeed, for small x, the Born 



xi. 



term is a function of xb^Q^. Hence, the choice /i^ = l/&^ is essentially equivalent to setting /i 

In ref. |2^, it was shown that the 6-space double logarithms exponentiate in higher orders. In the double 
logarithmic approximation, they give the suppression factor 

exp\-—CF\n'^(xQb)] (3.19) 

I An J 

for large b. The running of the coupling constant induces the next-to- leading logarithms (cf. p6| , p7t ). To get them, 
one should put as{k'^) = 47r/(/3olnfc^/A^) under the integral: 
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(^ p Ak±b± -1 

as Cf s{xQb) -^ — / cPk_L p as{k\ 



y^dy 



yy + kl/x^ 



(3.20) 



In general, the Sudakov effects are governed by the eikonal ||58| , p9| , p9[ (or cusp |6C|-|62f| ) anomalous dimension 

Cpa. 



-L cusp\Ois ) 



1 + ^ 

n 



,, , 67 7r2\ 5 ^^ 
^^(36-12-18^^ 



+ . 



(3.21) 



Clearly, only the as term of Tcusp{<3s) manifests itself in a one-loop calculation. To get further corrections p9| 



one should substitute Cpas/TT in Eq.(3.20) by Tcuspicts) and also use a two-loop expression for as{k'j_) and s{xQb) 



29|. Here, we restricted our analysis by the one- loop level. 



B. Vertex correction for the real photon 



For the real photon, the contribution of the vertex correction diagram 35 is given by 



Mb{x,Q,^,k^) = 



1 f ^ {x- OQ^ + xkl 



X fci [e(x - e)g2 + a;fc2 



e{i < x) . 



(3.22) 



Again, we concentrate on the term singular at k±^ = 0. It is convenient to split it into two parts. The first part 
is obtained by taking xQ^ from the {x — £,)Q^ term in the numerator and the second one by taking {—£,Q'^)- We 
represent the first part as 



fci K(a; - OQ^ + xkl 



eix 



ax 



(3.23) 



cQ^ + ki/u ax~o+xki/Q^ [^^Q2 + kimx-0Q^+xki] 



The last term here produces no divergences both for large and small k±. The 1/A:^ singularity is contained in the 
first term which we arranged to have a form of a product of the same Born term 1/(^Q^ -I- k\/a with a factor 
looking like a fc^-modified evolution kernel. Then we write this factor as a sum of a "plus" term and a 5{x — S) 
term: 



{e/x)B{^<x) 



i{x-0^xkl/Q^ \i{x-0+xkl/Q^ 



ie/x)e{^<x) 



+ 5{x-0 



{e/x)9{c<x) 

C{x-0+xkl/Q- 



dc. 



(3.24) 



As a result, the total contribution associated with the k± = singularity can be written as 



tr'ix,Q')^ d^ 



+ S{C - x) 



d^k^ j 1 f {e/x)9i^<x) 



S.Q^ + kll^ \ fci K^x - + xfci/Q2 



-5^{kA_) 



klUx-O+xkllQ^] 



d^k. 



kl 



C{x~C)+x~kl/Q^^ 



dC 



where the S'^{k±^) term comes from the second, "— ^Q^" part of the original expression ( 3.221 ). 
From this decomposition, we obtain the mass singularity term 



i oji < x) 

X (x-O J + 



L{bm) = Vb{£.,x)L{bm). 



the evolution-related contribution 



(3.25) 



(3.26) 
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Et{x,^;b) = - 






Ko{bQ^ax-C)/x) 



(3.27) 



and the Sudakov term 



Sb{x;bQ)^^ f S 



pik±b± _ -^ 



dC 



Cix ~C)+xkl/Q^^ 



= s{^/xQh) . 



(3.28) 



For large &, the latter behaves like 



Sh{x\bQ) 



■ In^ (V^Q6) 



(3.29) 



By analogy with Sa{x;bQ) which is a function of xQb we might expect that Sb{x;bQ) should be a function of 
xQb. Our calculation above shows that Sb{x;bQ) is a function of y/xQb. That this result is not unreasonable, 
can be justified in the following way. Note, that for small x, both the Born term B{x; bQ) and our Sb{x; bQ) are 
the functions of the same combination xb^Q^. Hence, integrating the product B{x; bQ)Sb{x] bQ) over b just gives 
1/Q^ multiplied by a constant factor: no Inx terms are produced. On the other hand, a In x term would appear 
if Sb{x; bQ) would behave like In {xQb) for large b. The explicit expression for diagram 35 given in ref. |q] has no 
In X terms. 



C. Box and self-energy diagrams 



In Feynman gauge, the box diagram 3c contribution in QCD 



Me(x,Q;e,fc±) = 



1 f x{^^Q' + kl) 



{x-C)^Q^ + xk 



5Q2 1 xki [eeo^ + ki] xki ux - 0Q2 + xki] 



0{i < x) 



(3.30) 



only by a numerical factor differs from that in a model with scalar or pseudoscalar gluons, in which Sudakov effects 
are absent. Hence, the k± — singularity produces only the evolution effects: 



1 



1 



M,ix;^,k^) = -^Vc{x,0 77^ + 



fci 



^Q^ 



(3.31) 



where Vc{x,^) is the relevant part 



Vc{x,O = -0{^<x) + ^d{^>x) 

X X 



(3.32) 



of the evolution kernel. Note, that Vc{x,£,) does not have a "plus" form by itself. The missing 6{x — £,) terms are 
provided by two quark self-energy diagrams 3d, e 



M, 



d+e 



1 1 /"^ 

xQ k^ Jo 



-9{c>x) + ^eic<x) 

X X 



dc = - 



1 



^X\a^ rv I 



S(x-0. 



(3.33) 



The third self-energy diagram 3/ has only the UV divergence: 



Mf 



1 



S{x-0 



C/x 



C{x-C)Q^ + xk 



2 d(^ d k±. 



(3.34) 



Combining evolution kernels from all the diagrams above, one obtains the total evolution kernel V{^, x) ( 2.23 ) 
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D. Standard vs. modified factorization 



Summarizing the findings of the previous subsections, we write the sum of the lowest-order term and one-loop 
diagrams in the impact parameter representation as 



F-y'jTroiQ ) — 



fli^^n^-hh^^-'^^ 



V{^,x) L{bm)+E{^,x;bQ) 



(3.35) 



+6{^ - x)S{x, bQ) + R{^, x; bQ) 



27r 



ip-nix) dx , 



where i?(^; bQ) is the 6- version of the Born term ( [3.13 ), V{S,, x) is the total evolution kernel, E{x, ^; bQ) is the sum 
of the evolution-related terms like ( 3.15| ), (3.27), S{x^ bQ) is the total Sudakov term given by Eqs.(3.17), ( |3.2§| ) and 
i?(f , X] bQ) accumulates all the remaining contributions coming from terms regular at fc_L = 0. Integrating over 
b and specifying the prescription for the renormalized distribution amplitude Lp.^{x]^), one would get the result 
(2.21) of the standard factorization scheme. In particular, the term i In x, most sensitive to the width of the 



distribution amplitude (/3^(x;/i), comes from a negative contribution — ^In x due to the Sudakov term S{x,bQ) 
and a positive contribution In x coming from the ?Tj-independent part of the convolution 



d^ I B{^-bQ)®V{^,x)L{bm)'^ 

J ^TT 



1 



+ Inx IniQ^/m^) -|- In^ x + f{x) 



(3.36) 



This convolution contains also terms denoted by f{x) which are less singular at x = 0. The total sum vanishes 
when integrated with the non-evolving asymptotic distribution amplitude ipT^(x). It does not vanish, however, when 
integrated with DA's differing from ip'^''{x). 

The logarithmic mass singularity Inm contained in the evolution term V{^,x) L{bm) is eliminated by absorbing 
it into the renormalized DA. The procedure used in the modified factorization approach of refs. | |29| , |30|] is to absorb 
ln(?nfe). As a result, one obtains the pion distribution amplitude (pnix; 1/6) normalized at the scale fi — 1/b. 
Making such a choice, one should realize that b is an integration variable and, to preserve the acquired precision, 
one must use the evolution equation to get (pTr{x; 1/6) for all relevant values of 6. In particular, if the distribution 
amplitude is assumed to have a CZ-type shape for large 6, it should be evolved towards the asymptotic shape for 



smaller 6 using Eq.(2.29). Modelhng (pTr{x; 1/6) by a function of x only amounts to neglecting the m-independent 
part of the convolution B{^]bQ) (g) V{£,,x)L{bm) ( 3.36| ). As noted before, this contribution contains In^ cc, hence, 
for extremely wide distribution amplitudes it can exceed that coming from the Sudakov term which only contains 
(-iln'x). 

In the formal 6 == limit, the function (/?^(a;;l/6) evolved according to the leading logarithm approximation 
formula (2.29), coincides with Lp'^^{x). However, the function E{x, ^; bQ) also develops a logarithmic singularity for 
small 6, because 

Ko {Qb...)^- ln(Q6) + . . . 



for small 6. Hence, two ln(6) singularities present in Eqs.(02),(3.35) compensate each other in the 6-^0 limit 
and the net coefficient in front of the evolution kernel is In^Q/m): the distribution amplitude evolves in fact only 
to the scale bmin ^ 1/Q corresponding to the resolving power of the external probe. Absorbing \n{Q/m) into the 
renormalized distribution amplitude one would get ip{x) —> ^{x; Q), with the large external momentum Q serving 
now as a factorization scale. Such a choice is usually made in the standard factorization approach, in which /j, is 
either a fixed constant, e.g., n — 1 GeV or proportional to the external momentum, ^ = aQ, with a being a fixed 
number. In particular, one can optimize the choice of the parameter a by taking the value producing the shape 
of ip{x; 1/6) averaged over the essential region of the 6-integration. Another point is that the pQCD evolution of 
ip{x] ii) is reliable only in a restricted region fj, }^ (Xq. Since the modified factorization involves integration over all 
6, we formally need to know the distribution amplitude f{x; 1/6) outside the perturbative region 6 ^ 1/mo- One 
should remember, however, that the Born term Ko{Qb^/xE) for finite x exponentially suppresses the large-6 region. 
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As a result, essential impact parameters b are '^ 1/Q- The suppression by the Born term disappears for small x 
when the effective scale becomes l/y^xQ^ rather than 1/Q. In this case, the suppression of the large-6 region is 
provided by the exponentiation of the Sudakov terms which is the crucial element of the modified factorization 
approach P9|,p0[ . As a result of the exponentiation, the series of [as In (Q6)]" terms, each of which tends to infinity 



as 6 ^ oo, is substituted by the exponential of Eq.( 3.19 ) type rapidly vanishing with growing b. Of course, for 



finite x, the Born term KQ{Qb\/xE) provides even stronger suppression of the large-6 region and the influence of 
the Sudakov factor is minor. Only for small x Sudakov effects become important. The relevant combination xQb 
in the Sudakov term of the diagram 3a converts into Qb, and the exponentiated Sudakov factor plays a primary 
role in squeezing the size of essential impact parameters. A special role of the small a;-valus in the &_L-integration is 
reflected by the — ^ hi x term resulting from the convolution of the Born term with the one-loop Sudakov factor: 

^ j B{x- Qb) S{x; Qb) <fb^ = -i^ ("-i In^ x - g{x) \ , (3.37) 

where g{x) stands for less singular terms. After integration with the asymptotic distribution amplitude, the 
(— i In X — g{x)) term gives approximately —9/4-1-0.05, to be compared with the magnitude —5/2 of the total one- 



loop correction (see discussion after Eq.( 2.27 )). Hence, the total one-loop correction in the case of the asymptotic 
DA is very close to the contribution of the Sudakov term alone (the deviation is only 12%). If the higher- loop 
corrections can be also approximated by the Sudakov contribution, then the exponentiated form would produce 
the all-order result in a rather compact form. 

Discussing the numerical significance of the Sudakov terms, we should keep in mind that all the logarithmic en- 
hancements In [Qb) are perfectly integrable and that the region of small x, where the Sudakov terms are important, 
is small itself: after 6x- and x-integrations, there are no especially large contributions in the final result. The total 
one-loop correction is only about 20%. Hence, the exponentiation of the Sudakov terms would alter the one-loop 
corrected result for the form factor by just a few percent, which is similar to the accuracy of approximating the 
total contribution by the Sudakov term at one loop. Note also that a few percent change may be smaller than 
the contribution generated by the one-loop terms E{x,£,,Qb), R{x,£,,Qb) and the effects due to the 6-dependence 
of the renormalized distribution amplitude ipTr{x] 1/6). Moreover, for a wide DA, the latter are comparable to or 
exceeding the Sudakov contributions. In principle, one can try to explicitly include these corrections within the 
MFA framework, but the result would not have a simple form anymore. In this situation, instead of dealing with 



convolutions of Bessel functions, one may prefer to use the result (2.21) of the standard factorization approach 
which has a simple form with easily controllable accuracy. Another bonus of using the SFA is the ability of (pTr{x; Q) 
to fully absorb the necessary nonperturbative information: increasing Q we do not need to make any assumptions 
about the shape of (/^^(a;; ^J-) at smaller values fi < Q oi the factorization scale fi. 

IV. INCLUSION OF PRIMORDIAL TRANSVERSE MOMENTUM 

A. Brodsky-Lepage interpolation 

Despite our persistent efforts, we failed so far to find any traces of contributions capable of producing a series of 
transverse-momentum-related power corrections to the leading pQCD result. Recall that we investigated first the 
higher-twist contributions due to operators with contracted covariant derivatives D'^ . . . D^ which are the standard 
candidates to describe the /cj^-effects in the OPE-like factorization approaches. We observed that, for the simplest 
handbag diagram, these operators do not produce the expected infinite chain of (1/Q^)" power corrections. Then 
we studied one-loop radiative corrections in the Sudakov and impact-parameter representations. Our results are 
in full accord with the corresponding expressions of the MFA p9| , p0|| . But they also completely agree with the 
one-loop results 0-0] of the SFA, i.e., they do not contain any power corrections. Nevertheless, F^.^7r<'(Q^) ^ VQ^ 
cannot be a true behavior of F^-^^o^Q"^) in the low-Q^ region, especially since the Q^ = limit of FY-yaoiQ"^) is 
known to be finite and normalized by the 7r° -^ 77 decay rate. The value of F^-.^^o{Q) in QCD Eo] is fixed by the 
axial anomaly 
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^7*7-°(0) = r^- (4.1) 



If the shape of the pion DA is specified, the large-Q^ behavior is also known. For the asymptotic DA, 

J^,Vo(Q2) = ^. (4.2) 

Long ago, Brodsky and Lepage H] proposed the interpolation formula 



which reproduces both the Q^ = value (O) and the high-Q^ behaviour given by Eq. (12). The BL-interpolation 
formula ( |4.3| ) has a monopole form with the scale sq = 47r^/J « 0.67 GeV^ numerically close to the p-meson mass 
squared: nil w 0.6 GeV^. Thus, the BL-interpolation suggests a form similar to that based on the VMD expectation 
Fj»y^o{Q'^) = l/[7r/7r(l + (5^/mp)]. In the VMD-approach, the p-meson mass Wp serves as a parameter determining 
the pion charge radius, and it is only natural to expect that the tower of (so/Q^)^-corrections suggested by the 
BL-interpolation formula can be explained by intrinsic transverse momentum effects. The only problem is how 



to get Eq.(4.3) (or anything similar to it) from QCD, i.e., how to construct an expression which would provide a 
good model both in perturbative and nonperturbative regimes. Before proposing our variant of the solution to this 
problem, let us discuss briefly two recent attempts |2l| , |26| to include intrinsic transverse momentum effects into 
the description of the 7*77r'^ form factor. 

B. Extrapolation of perturbative results 

As emphasized above, despite the fact that the denominator of the Born term 1/(^(5^ -t- k±/£,) is fcj_-modified 
compared to its collinear approximation £,Q^, convoluting B{^; bQ) with S{£,; bQ) one would enjoy no power modi- 
fications of the canonical 1/Q^-behavior, i.e., the transverse-momentum effects included in the Sudakov term and 
other one-loop corrections do not correspond to any higher- twist contributions. The obvious reason is that, apart 
from the IR regulator mass m (producing a logarithmic dependence Inm which is absorbed into ^{x; /x)), the large 
momentum Q is the only scale that appears in the relevant fc^-integrals. 

In general, the fact that some contribution is written as an integral over the transverse momentum kj_ or the 
impact parameter b± does not necessarily mean that something beyond the leading twist is included. To illustrate 



this point, we note that even the lowest-order, "purely collinear" contribution (2.11) can be written in the impact 



parameter representation. A possible form is suggested by the one-loop calculation: 



FoiQ^) = \ I dx f xKo (V^x62Q2^ ip^ix)dH, (4.4) 



where xKq{^ xxb'^Q'^) is the impact-parameter profile of the modified propagator l/{xQ^ + k\/x) (see Eq.(^.13)) 



Though the 6- version of the quark propagator explicitly depends on b, integrating over b in Eq.(L4) gives a simple 
power result l/Q^ without any subleading power corrections. This phenomenon can be traced to the absence of 



the 6-dependence in the distribution amplitude. In the momentum representation, Eq.(L4) is equivalent to using 
tpT;(x)5'^{k\^ for the qqn vertex: 

^2^ „ 47r [^ f ip^{x)5'^{k^_) ^ 



However, as we have seen in the preceding section, radiative corrections generate terms with less trivial fcj^- 
dependence. In particular, the one-loop correction contains a^/fc^ terms. As a result, the fc^-dependence of the 
qqir vertex at one loop is 
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^-(OS\k^) + j^^ J^ Vi^,x)ipAx)dx + .... (4.6) 

In the impact parameter representation, the sum of S^{k±) and 1/fc^ terms is converted into a more suggestive 
combination 

^AO - 7^ Hbm) / V{^, x) ^4x) dx, (4.7) 

which can be understood as the two first terms of the a^-expansion of the expression for the leading-logarithm 
evolved distribution amplitude tp{^, \/b) written symbolically as 



exp 



\n(hm)V 

27r ^ . 



® V5. 



Since all the conclusions made from the studies of one-loop corrections are based on perturbative analysis, strictly 
speaking, they are only applicable to transverse momenta which are large enough fj. Furthermore, there are no 
special reasons to expect that formulas derived for momenta k^ generated by perturbative gluon radiation are still 
true in the small-fcj^ region dominated by primordial (or intrinsic) transverse momentum. Still, it is tempting to 
extend the leading-logarithm convolution formula 

F(Q2) = I / dx I KoiV^Qb) ip{x; 1/6) d% (4.8) 

into the nonperturbative region. To do this, we should substitute the distribution amplitude 'p{x; 1/6) by a function 
which reflects (or models) the nonperturbative 6-dependence. 

In the light-cone approach Q, the basic object is the wave function \['(x, fcj_) which depends both on the fraction 
variable x and transverse momentum k±. In QCD, it is customary to split ^(x, k±) into two components. The soft 
component ^*°'^*(x, A:_l) is due to the nonperturbative part of the QCD interaction and its width is determined 
by the size of the relevant qq bound state. It is expected that ^^"■'^'(a;, fcj^) rapidly {e.g., exponentially) decreases 
for large fcj^. In our perturbative lowest-twist treatment above, the soft wave function '^''°f*{x,k±) was imitated 
by (pTT{x)S'^{k±). The pQCD interaction (gluon radiative corrections) produces the hard component ^'"""^(x, /cj_) 
which behaves like as/fc^ at large kj_. The distribution amplitude (pTr{x) can be treated as the integral of the wave 
function $(a;,fcx) over k± (cf. [Q): 

VA^)^^ j^{x,k^)d^k^. (4.9) 

For '^^°^*{x,k±), this integral perfectly converges. However, the perturbative l/k\-iai\ generates logarithmic 
divergences. Hence, one should supplement this definition by some regularization procedure specified by a cut-off 
parameter /j,: tfTr{x) -^ (pTr{x, fi). The "cut-off" should be understood in a broad sense. It may be imposed literally 
k\ < fp or one can use more gentle procedures based, say, on dimensional regularization. In other words, (frrix) 
is a scheme-dependent object: (/J,r(x) — *■ (/jV (x). The choice of a specific scheme S' is a matter of convenience. In 
particular, the Fourier transform 

$(x,6) ^ ^— j e-'^^''^-^{x,k^)d^k^ (4.10) 

to the impact parameter representation can also be treated (at least, for small 6)rj as a regularization scheme for 
the integral defining the distribution amplitude: 



^^In particular, speaking about the double-logarithmic behavior "at large h" we imply that 6 may be much larger numerically 
than l/Q but is still within the pQCD applicability range. 

'■■'■ The basic difference between tp^{x\ 1/6) and *I'(a;, h) is that f ^ti{x\ /i) is given by the same constant /^ for any /x while 
r '^{x, h)dx in general depends on b. 
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^i^\x■^i^l/h)^^^{x,b) ■ b^O. (4.11) 

This observation suggests the extrapolation of the convolution formula into the nonpcrturbative region by sub- 



stituting tp{x;l/b) in Eq.(4.8) by the 6-space wave function '^{x,b) (see ref. [^). Since the fcj^-effects are only 
essential when xQ^ (i.e., x) is small, one can either use the original combination \/xEQb in the argument of the 



Born term Ko{^/xxQb) or substitute it by yJxQb. In particular, a modified version of the convolution formula (4.8) 
written in the fcj_-representation 

^ '■ *(a;,fci^ 



F,.,..m--^^li^l^§^/^.. (4.12) 

is the starting point of the analysis by Jakob et al. pi[ . In this expression, a simpler form xQ"^ + k\ is used for 
the modified denominator of the "hard" quark propagator instead of the combination xQ^ + k\/x which appears 



in our Eq.(3.13). However, since the difference is proportional to k\ and vanishes for a; = 0, the two forms have 
essentially the same footing. As a model for "^{x^k^^), Jakob et al. pl| use the ansatz [p| with the exponential 
dependence on the combination k\lxx (or Gaussian dependence on fcj^). We write it in a form similar to that used 
in ref. ^■. 



crV6 XX \ 



y\ 



2axx J 



(4.13) 



where a is the width parameter and fTr{x) is the desired pion distribution amplitude^ In the 6j^-representation, 
the model wave function is 

¥^\x,b^) = ^(^^(x)exp (~]^bi(Txx\ . (4.14) 

The model is restricted by two conditions taken from ref. 1^]. First, the two-body Fock component of the pion 
light-cone wave function ^(a;, k±) is required to satisfy the constraint 

>"=/*(-.'=jS^ = ^ (4.15) 

imposed by the tt ^^ ^u rate. This gives the usual normalization condition for the pion DA 



ip-^{x)dx^U. (4.16) 







The second condition specifies the value of the x-integral of ^(x, k±) at zero transverse momentum 

1 /^ 

^{x,k±^Q)dx ^^. (4.17) 

J-rr 



For the model ansatz (4.13), this condition results in the following constraint for the /o-integral 



/o = — / ^^{x)—^—. (4.18) 

Jtt Jo X So 



^^In the original model pi k']_ appears in the combination k']_ + M^ where Mq is the constituent quark mass. As a result, 
the distribution amplitude (pTr{x) is exponentially suppressed like exp[—Mq /2axx\ in the end-point regions. Jakob et al., 
however, follow Chibisov and Zhitnitsky |63] who insist that the constituent quark mass Mq should not appear in QCD- 
motivated models for '^{x,k±_). In particular, Mq does not appear in the model wave function '^^^^'{x,k±_) P3] based on 
local quark- hadron duality (see Section V below): only the current quark masses m, (usually set to zero for u and d quarks) 
are present in QCD Feynman integrals. 
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In obtaining Eq.(4.18), we incorporated the symmetry property ipTr{x 
the notation sq for the important combination 47r^/^ 



Since 11^" 



3 and I^^ 



iPty{x) of the pion DA and used again 
5, the width parameters are 
cr"'' ^ So- 0.67 GeV^ and a^^ = |so « 1.11 GeF^. 

In the form (4.17), the second condition was derived in ref. |2] from the requirement that the tt" -^ 77 decay 
rate (or, what is the same, F^.jt^o{Q^ = 0) ) calculated within the light-cone approach coincides with that given by 
the axial anomaly. It is easy to see, however, that in the Q^ ^ limit, the fc^-integral in Eq.(4.12) logarithmically 
diverges in the small-fcj_ region for any function which is nonvanishing at kj_ = 0. Note, that 'I'(x, /cj^ = 0) 



cannot vanish if we wish to satisfy the condition (4.17). Rather ironically, the condition which presumably should 
secure the correct value for Fj*yT^o{Q^) at Q^ = guarantees instead that the extrapolation formula diverges at 
that point. This gives a clear warning that one should be very careful using the simplest extrapolation: it is 
difficult to judge a priori how reliably the formula failing for Q^ = models the subasymptotic effects for moderate 
Q^. The authors of ref. pl| ] also include the Sudakov exponential in which they take a symmetric combination 
s{xQb) + s{xQb). As noted earlier, our one-loop calculation in Sect.IIIi? shows that for F^..^^o((5^) one should use 
s{y/xQb) instead of s{xQb). Our final observation is that expanding Eq.(4.12) in k'^/Q'^ one would get an infinite 
series of power corrections under the x-integral. According to our general result, the handbag diagram should not 
produce a chain of higher-twist contributions. Hence, the extrapolation formula cannot be interpreted simply as a 
transverse-momentum-corrected expression for the handbag diagram. 



C. Transverse momentum in the light-cone formalism 



Another attempt to model the subasymptotic corrections was made in ref. Ig^l • It is based on the Brodsky-Lepage 
formula lit] for the two-body (i.e., qq) contribution to the 7*777*^ form factor in the light-cone formalism: 



{e^xq^)F^^JQ') 



tt^VG Jo 



dx 



(xq± + k^y — Id 



(4.19) 



Here, g_L is a two-dimensional vector in the transverse plane satisfying g^ ~ Q^ , e± is a vector orthogonal to q± and 
also lying in the transverse plane M and the cross denotes the vector product. Again, the wave function is chosen 
in the Gaussian form ( [1.13| ) satisfying the constraints ( 4.15| ) and ( 4.17 ) *** . Though the integrand of Eq.( 4.19| ) 
looks rather singular, there are no problems with the convergence of the fc_L-integral in the q± ^ limit. The 
result is finite, since 



ql+kl 



(qa 



k,V-' 



Tr6^{kA 



(4.20) 



gj_-*0 



for any test function ^(x, fcj^) which depends on fcj^ through k^. Because of the S (kj 
is determined by the wave function at zero transverse momentum. 

In ref. p6|, it is claimed that the k±/Q expansion of Eq.(4.1£) produces large 
the leading- twist result. In fact, when ^{x,k±) has an exponential fc^-dependence, it is trivial to calculate the 
fc^-integral explicitly 



)-function, the Q^ = result 
■'higher-twist" corrections to 



F^^^^oiQ^) - Y / 



xQ^ 



1 — exp 



2xa 



dx 



(4.21) 



***As emphasized recently by KroU [Bl|, Cao et al. use constituent quark masses M, 
exponential suppression exp 

their "CZ" model is 3.71 rather than 5, i.e., despite zero at a; = 1/2, such a model gives a rather narrow DA, which is closer 
in this sense to the asymptotic DA rather than to the original CZ one. 



330 MeV which produces a strong 
Mq/2axx] of the end-point regions. As a result, the /-integral for the DA corresponding to 
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to see that the correction term in the integrand of Eq.(4.21) has an exponentiaUy decreasing rather than a power 
behavior for large Q^ . This result agrees with our general statement that the handbag diagram contains no higher- 
twist contributions. Our analysis works in this case since the Brodsky-Lepage formula ( 4.19| ) corresponds to the 
handbag contribution written in the light-cone variables without any approximation. Just like in the covariant 
treatment, the naively expected series of power corrections ((fcj^)/Q^)" does not appear because the expansion of 



xq\^ -I- /c^ 



(4.22) 



contains only traceless combinations. Indeed, multiplying ( 4.22 ) by q^/Q"^ and defining (fcj^g^) = jfej^jQcos^, we 
obtain 



xQ"^ + |A;^|Qcos( 



1 



QV x'^Q'^ + 2x\ki_\Q cos (l) + k\ 



l^\e{\k^\<xQ) + Y,{-ir 



(4.23) 



M\\{\k^\<xQ)-(^^] e{\k^\>xQ) 



cos(n0) 



For a wave function ^(x, fcj^) depending on k± through k\ only, all the oscillating terms proportional to cos(n(/>) 
(i.e., to Chebyshev polynomials T„(cos(/i) corresponding to traceless combinations in two dimensions) vanish after 
the angular integration. Only the n = term written outside the sum over n gives a nonzero result. Hence, for 
the wave functions of ^{x, k±) — ^l^ix, fc^) type, we can write 



-f'~.^jro(V ) 



dx 



nVe Jo xQ 



xQ 



ip{x, k^) k±dkj 



(4.24) 



This means that the leading 1/xQ^ term in Eq. ( [4.21D com es from the integral over all fc^'s while the exponential 
correction appears because the integration region in ( [4.24 ) is restricted by k± < xQ. Another subtlety is that the 
Q2 _ value 



F'3'i{Q^ ^ 0) 



2TTf^ 



dictated by eqs.(4.I7) and (4.20) (and manifest in Eq.(4.21) ) gives only a half of what is needed to get the correct 
7^0 _s. ry^ j.a,te (4.1) . As explained in ref. |2|, the other half comes from the term which can be interpreted as 
the contribution of the qq^ Fock component of the pion wave function. In a formal pQCD diagrammatics, this 
contribution is represented by graphs containing the gluons coupling to the quark line between the photon vertices. 
For high Q^, such diagrams correspond to higher- twist corrections associated with the qG . . . Gq operators. In this 
sense, the result of ref. is equivalent to a nonperturbative constraint on the Q^ ^ limit of such contributions. 
One can expect that the qqj contribution decreases as 1/Q^ or faster for large Q^ since it contains higher twists 
only. Interpretation of this contribution in terms of the qqj Fock component is restricted to the case of real 7: 
ref. H gives no expression beyond the Q^ — point. In ref. pq] this contribution is not included. However, if the 
terms which double the result for Q^ = are not included, it is premature to make specific quantitative statements 
about the size of subasymptotic corrections in the region of moderate Q^. 

We may also wonder why the formulas ( 4.12| ) and ( 4.19 ) corresponding to two attempts to include the primordial 
transverse momentum have such a strikingly different analytic structure. In particular, the denominator of the 
integrand of Eq.(4.19|) vanishes for fcj_ = —xq± while that of Eq.(4.12) is finite for all k±_ provided that qj_ ^ 0. 



The answer is very simple: the two expressions imply two different definitions of what is longitudinal and what 
is transverse. Eq. (4.12) is based on the Sudakov decomposition in which the momentum qi of the real photon 
has only the light-cone "plus" component while the momentum p of the pion has only the light-cone "minus" 
component. As a result, the momentum transfer (72 — P^ li in the Sudakov variables is purely longitudinal and has 
both plus and minus components, with g| = —2{qip). On the other hand, the Brodsky-Lepage formula corresponds 
to the infinite momentum frame in which the plus components of gi and p coincide. The plus component of the 
momentum transfer (72 vanishes in this frame, but (72 has a nonzero transverse component q±, with \qj_\ = Q or 
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^2 = ~Q±- Evidently, the two frames cannot be obtained from one another by a boost. Furthermore, one should 
not expect a diagram by diagram correspondence between the two approaches. The main purpose of imposing the 
requirement (73 = in the light-cone approach is to avoid the Z-graphs. However, in Sudakov variables (and in any 
approach in which 92 has a non-zero plus component) the Z-graphs should be added to reproduce the light-cone 
result (cf. 11). 

Both the approaches pi|j26[ discussed above fail to reproduce the Q^ = value corresponding to the axial 
anomaly. Our point of view is that complying with the anomaly constraint should be a minimal requirement for 
any model of subasymptotic effects in the 7*771" form factor. A maximalist attitude is that such a fundamental 
constraint should be satisfied automatically rather than imposed as an external condition. This can be only realized 
in an approach which is directly related to QCD and produces anomaly as a consequence of QCD dynamics. 

V. QUARK-HADRON DUALITY AND EFFECTIVE WAVE FUNCTION 

A. QCD sum rule calculation of /,r and local duality 

QCD sum rules provide us with the approach which deals both with perturbative and nonperturbative aspects 
of QCD. The basic idea of the QCD sum rule approach [p5| is the quark-hadron duality, i.e., the possibility to 
describe one and the same object in terms of either quark/gluon or hadronic fields. To get information about 
the pion, the QCD sum rule practitioners usually analyze correlators involving the axial current. In particular, to 
calculate /^ one should consider the p^p^-part of the correlator of two axial currents: 

n^"^(rt = I J e'P-{0\T (J5^(a;) J5.(0) )| 0) d^x = p^^p^n^ip^) - g,,,ni{p'). (5.1) 

The dispersion relation 

n2(p^) = -/ \, ds+ ''subtractions'' (5.2) 

represents Il2{p^) as an integral over hadronic spectrum with the spectral density phadron^g-j determined by pro- 
jections 

(0|j5;.(0)k;P)-*/.P^, (5.3) 

etc. , of the axial current onto hadronic states 

phadron^^^ = TTf^S{s - m^) -|- 7r/^/(s - m\j + "higher states" (5.4) 

( f^^^ ~ 130.7 A/eV in our normalization). On the other hand, when the probing virtuality is negative and large, 
one can use the operator product expansion 



n,{p^) ^ nr^'^ip^) + ^{asGG) + ^asimf + ... (5.5) 



where II'^""'^ {p^) is the perturbative version of 1X2(15^) given by a sum of pQCD Feynman diagrams while the 
condensate terms (GG), (qq), etc. (with perturbatively calculable coefRcients A,B, see Eq.( ^.10D below), de- 
scribe/parametrize the nontrivial structure of the QCD vacuum. For the quark amplitude XIj""'^ (p^); one can 



also write down the dispersion relation (|5.2D , with p{s) substituted by its perturbative analog p'""'' (s): 



P^""'''(5)-^(l + ^ + ...) (5.6) 



(we neglect light quark masses). Hence, for large —p, one can write 
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1 r^^ Jxadronf „\ „quark / „\ A r? 

i / '- ^^^"f ^'^ ds . 4(a.GG) + ^^aM,? + ■■■■ (5.7) 

This expression essentially states that the condensate terms describe the difference between the quark and hadron 
spectra. At this point, using the known values of the condensates, one can try to construct a model for the hadronic 
spectrum. In the axial-current channel, one has an infinitely narrow pion peak p^ = Trf^d{s — rn^), a rather wide 
peak at s w 1.6 GeV^ corresponding to Ai and then a "continuum" at higher energies. The simplest approximation 
is to treat Ai also as a part of the continuum, i.e., to use the model 

phadron^^^ ~ 7:^8 - ml) + p'^^^'''^ {s) 0(s > Sq), (5.8) 

in which all the higher resonances including the Ai are approximated by the quark spectral density starting at 
some effective threshold Sq. Neglecting the pion mass and requiring the best agreement between the two sides of 
the resulting sum rule 

/•2 1 /-So „quark( \ A f> 

J^^_ P ^^^, ^ -as{GG) + ^a,{qqf + . . . (5.9) 

P ^ Jo s — p^ p^ p° 

in the region of large p^, we can fit the remaining parameters /^ and sq characterizing the model spectrum. In 
practice, the more convenient SVZ-borelized version pq] of this sum rule 

is used for actual fitting. Using the standard values for the condensates {GG), {qq}^, the scale Sq is adjusted to 



get an (almost) constant result for the rhs of Eq.(5.10) starting with the minimal possible value of the SVZ-Borel 
parameter AP. The magnitude of /,r extracted in this way, is very close to its experimental value /^^^ « 130 MeV. 
Of course, changing the values of the condensates, one would get the best stability for a different magnitude 
of the effective threshold sq, and the resulting value of /^r would also change. There exists an evident correlation 
between the values of /,r and sq since, in the Af ^ — > oo limit, the sum rule reduces to the local duality relation 

f^ = - r p''^-^Hs)ds. (5.11) 



Using the explicit lowest-order expression Pq""'^ (s) ~ l/47r, we get 

so = 47rV'. (5.12) 

Note that sq — 47r^/^ coincides with the combination which appears in the Brodsky-Lepage interpolation formula 



(4.3) 



B. Quark-hadron duality for the F^,^,^o \Q) form factor 

Information about the 7*7* — > tt" form factor can be extracted from the three-point correlation function ||66|1 

.^0^.(91,92) ^^j d^xd^y e-^«^^-^«^«(0|T{J^(x) J,(2/)j5o(0)}|0) (5.13) 



calculated in the region where all the virtualities q\ = ~q^, g| = —Q^ and p^ = {qi + (72)^ are spacelike. 

The form factor F^.,^t^o{q^ ^Q"^) appears in the invariant amplitude Fyp^^q^^Q^j corresponding to the tensor 
structure e^iypcrPaqiq2- The dispersion relation for the three-point amplitude 

F(p^,q^,Q^) ^ - '- \; ' ds+ ''subtractions'' (5.14) 

-^ Jo s-p^ 
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specifies the relevant spectral density p [s, q^ , Q^)- For the hadronic spectrum we assume again the "first resonance 
plus perturbative continuum" ansatz 

phadron (^^ ^2^ g2) ^ nUF^,^,^o{q^ Q^)5{s - m^ ) + 0(s > Sq) p'^^'^^'^s , q\ Q^) . (5.15) 

The lowest-order perturbative spectral density p'^'^°-^^{s, q^, Q^) is given by the Feynman parameter representation 

P«"'^'^'=(S, q^,Q^) =2 Tsfs- q^XiX3 + Q^X2X3 \ s(l_J2 xA dxidX2dX3 . (5.16) 

Scaling the integration variables: xi +2:2 = y, X2 = xy, xi = (1 — x)y = xy and taking trivial integrals over x^ and 
y, we get 

p«-'=(., q\Q^) = 2 /' ^!(;^' +f !)' dx . (5.17) 

Jq [sxx + xQ-^ + a;<7^J-^ 

The variable x here can be treated as the light-cone fraction of the pion momentum p carried by one of the quarks. 



In particular, the denominator of the integrand in Eq.(5.17) is related to that of the hard quark propagator: 
{qi — xpY ~ —{xQ'^ + xq^ + sxx). 

Putting one photon on shell, g^ = 0, we can easily calculate the x-integral: 

This result explicitly shows that if the larger virtuality Q^ also tends to zero, the spectral density p''"°'''"'(s, Q^) 
becomes narrower and higher, approaching 5{s) in the Q^ ^ limit (cf. p7[|). Thus, the perturbative triangle 
diagram dictates that two real photons can produce only a single massless pseudoscalar state: there are no other 
states in the spectrum of final hadrons (cf. pq]). As Q^ increases, the spectral function broadens, i.e., higher states 
can also be produced. 

A detailed study of the QCD sum rule for the Ky.^^o(Q^) form factor was performed in refs. [g^,^. The results 
of this investigation are rather close to those based on the simple local quark-hadron duality ansatz: 



■^U Jo 



(5.19) 



Using the explicit expression for p*"'"''''(s, Q^), we can write 

■^U Jo Jo [sxx + xQ'^]^ "" -n-fTri^ + Q^/so) 



F^h^JQ') = ^Cdx r ,„!y^2L ds = ,, ,, ' ,,^, . (5.20) 



This result coincides with the Brodsky-Lepage interpolation formula (4.3) 



C. Effective wave function 

The formulas based on the local quark-hadron duality prescription can be interpreted in terms of the effective 
two-body light-cone wave function 123]. Consider the lowest-order perturbative spectral density for the two-point 
correlator. It can be written as the Cutkosky-cut quark loop integral 

^,na./c(^) ^ _L I ^ (^1 _ ^^ e[k+) 5 (fc2) e{p+ - k+) 5 {{p - kf) d'k (5.21) 

where s ^ p^. Introducing the light-cone variables for p and k: 

P^{p+ = P,p- = s/P,p^ - 0} ; k^{k+= xP, k^,k^} 
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and integrating over fc„ , we get 



.--(.) = if ^./^(.s- I )^^fc.. (5.22) 



The delta- function here expresses the fact that the Ught-cone combination k^/xx coincides with s = p^, the 
invariant mass of the qq pair. Substituting this expression for p''"'"''^(s) into the local duality formula (5.11), we 
obtain the following representation for /^ 

fl = ^ I dx I e{kl< xxs^) <fk^ . (5.23) 



It has the structure similar to the expression for f.^ in the light-cone formalism [|l| (cf. ( 4.15 )) 



U^V6 J dxj^ix,k^)^. (5.24) 
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To cast the local duality result ( 5.23| ) into the form of Eg. ( ^.24 ), we introduce the "local duality" wave function 
for the pion: 

*^^(x, fci) = ^ e{kl < xxso) . (5.25) 

In 

The specific form dictated by the local duality implies that ^^^(x, k±) simply imposes a sharp cut-off at k'j_xx = sq. 
In the 5^-space, the effective wave function can be written as 

*^^(a;,5i) = ^-— V^I5^Ji(&±\/Ssi^), (5.26) 

where Ji{z) is the Bessel function. 



D. Effective wave function and F^,^^a{Q^) form factor 



Consider now the local duality expression (5.20|) for F^-,^^o{Q'^). Replacing s, the invariant mass of the qq pair, 



by its light-cone equivalent k\/xx^ we get F^P o(Q^) as an integral over the longitudinal momentum fraction x 
and the transverse momentum fc^: 

F^,%.m -^J^-J (^|?Tw '^'" - ""'°^ "'"" • ^'-'"^ 



Now, introducing the effective wave function '9^^{x,k±) given by ( ^.25 ), we write F^^ {Q^) i^ the "light-cone 
form" : 



^"-"'«'> -j^-jy- / ijiTw*'°(-'j^^'- 



(5.28) 



In the impact parameter representation, this formula looks like 

1 
F^P^^oiQ^) = -^ f dx f xQ^b'' K2 {V^bQ) *^^(x, 6^)^26^. (5.29) 



The function K2 (^/xbQ), where K2{z) is the modified Bessel function, originates from the new version of the Born 
term written in the 6-space 
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B{x; bQ)^^l e-^-^- J^^^^ '^"k^ = \ -Q"b"K^ i^^Q) ■ (5-30) 

Note that B{x;bQ) is finite for 6 = 0: B{x;0) = 1 while the "old" Bom term B{x;bQ) = xKq (y/^bQ) ( |3.13| ) 
has a logarithmic singularity at the origin of the &-space. The expression ( 5.28 ) looks similar to the extrapolation 
formula (4.12|). Furthermore, since 



W? _ 1 2fci , fcl ^^^^^ 



2 A3 



{xQ^ + kl)3 xQ^ + kl (xQ^ + klf (xQ^ + kl) 

the two fc_L-modifications of the hard quark propagator 1/xQ^ differ only by O(fc^) terms invisible in the analysis of 
effects induced by the 1/A;^ singularity at small k±. However, this difference is very essential when one extrapolates 



into the region of small Q^. To demonstrate this, let us analyze Eg. ( 5. 28 ) in some particular limits. For real photons, 
using the fact that 



M^ 1 



-6{ki) (5.32) 



(/x2 + klf 2 

in the /i^ — > limit, we obtain that the tt^ -^ 77 decay rate is determined by the magnitude of the LD wave 
function at zero transverse momentum: 

^r7-°(0) = ^/ *^''(a^,^x=0)dx. (5.33) 

This requirement is similar to that in the Brodsky-Lepage formalism. However, according to the explicit form 



(5.25) of '^^^{x,kj_ = 0), the integral ( [3.33| ) is twice larger than the constraint (4.17) imposed on the valence 
qq light-cone wave function. As a result, the local duality formula exactly reproduces the F^,^t^o[0) value (4.1) 
dictated by the axial anomaly. This outcome can be interpreted by saying that 'i'^^{x,k±) is an effective wave 
function (cf. |69| ) describing the soft content of all qG . . . Gq Fock components of the usual light-cone approach 
(see also [|3|). Note, that higher-order radiative corrections to the perturbative spectral density p'^^"-'^^ {s , Q"^) are 
explicitly accompanied by the as(/i|j)/7r factors per each extra loop. After integration over the duality interval 
< s < So, there are two physical scales: sq and Q^ . At low Q^, the duality interval Sq sets the scale at the low- 
momentum end of the f/T^-divergent integrals, hence, a natural choice for the normalization scale ^n is /i|j ~ sq. 
At high Q^, the short-distance dominated parts of the higher-order corrections should reproduce the pQCD results 
which suggest /i|j ^ Q^ for these terms. In any case, suppression by at least as{sQ)/TT ^ 0.1 per each extra loop is 
guaranteed. Since sq 3> A^, the gluonic corrections to p'^'"""'^(s, Q^) are suppressed by powers of as{sQ)/'iT ^ 0.1. In 
other words, the higher-order diagrams contributing to p'"'^'''^(s, Q^) correspond to exchange of hard gluons whose 
wave lengths are larger than 1/ ^/sq. 

When Q^ is so large that the fc^-term can be neglected, we get the expression 



F^,%om ^'^j ^^D{x,k^)<fk^ +0(1/0^). (5.34) 



"^ dx 

Identifying the wave function integrated over the transverse momentum with the pion distribution amplitude 

f^'^ix) = -^ J vl/^^(a:, k^) d?k^ = 6/,x(l - x), (5.35) 

we obtain the lowest-order pQCD formula ( ^.ll| ) 

F^>^AQ^)W-^o. = ^ f^^dx + 0{l/Q') (5.36) 

for the large-Q^ behavior of the 7*7 -^ 7r° transition form factor. 
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To summarize, the local duality formula ( 5.20| ) exactly reproduces the Brodsky-Lepage interpolation (4.3) be- 
tween the Q^ = value I/tt/tt fixed by the axial anomaly and the leading large-Q^ term 4:Trf^/Q'^ calculated for 
the asymptotic form of the pion distribution amplitude. 

The application of the local duality ansatz in a general situation when both photons are virtual was discussed 
in ref . P3[ . The basic formula written in terms of the effective wave function is given by 



pLD ( „2 ^2 



1 



--"'•« = ^/. ''"•"(■■'"•«^» = ^/.^7 K+.Q» + .,^P 



xxixQ"^ + xq^Y .T.LD 



*^^(a;,fc_L)d A:j 



(5.37) 



For g = Q = it satisfies the anomaly constraint (4.1), while when both q and Q are large it reduces to the 
pQCD formula 



Fj'f'TT°{Q )lg^Q2^oo 



4n 



xQ^ + xq^ 



dx + 0{\/Q^). 



(5.38) 



E. Extended local duality 



Note, that the pion distribution amplitude ( 5.35 ) produced by the local duality prescription coincides with the 
asymptotic DA. To model wave functions corresponding to DA's different from (/?"''(a;), we propose to use the sharp 
cut-off analog of the Gaussian model ( 4.13| ) : 



* 



(LD) 



(x,k±_) 



a\f% XX 



(fcj^ < xxa) 



(5.39) 



where a is again the width parameter and ^-nix) the desired DA, which satisfies the standard /^-normalization 
constraint ( 4.16 ). To guarantee the anomaly result for the tt^ — > 77 rate, we impose the following constraint on 
the a;-integral of '^^^^'(x, k±) at zero transverse momentum 

/' vl,(ic)(a:,fc^=0)dx =^. (5.40) 

Jo Jt 

Substituting the model ansatz ( ^.39| ), we derive from this constraint the condition for the standard integral /q 

dx 3(7 



u 



iPrr{x)- 



SO 



(5.41) 



where Sq is the basic combination sq — in'^f^. Taking /q* — 3 and Iq^ = 5, we fix the width parameters cr°^ 



So 



and 



^cz 



fSO 



l.llGeV^. Note, that in the CZ calculation [|3|, the duality interval was 0.75 GeF^ for the 



zeroth moment of the DA and 1.5 GeV"^ for the second one; our effective duality interval a'^^ for the CZ-type DA 
appears to be the average of these two. Using the ansatz ( 5.39| ) in Eq.( ^.2g ) and integrating over the transverse 
momentum, we obtain 



F^.^^o(Q^) - ~^ 



2vr [^ ifT^ix) 



1- 



1 



{l + xa/Q^y 



dx. 



(5.42) 



This formula has correct limits both for Q^ — and large Q^. For the asymptotic distribution amplitude. Eg. ( 5. 42 ) 
produces the expression (5.20) coinciding with the Brodsky-Lepage interpolation formula. For the Chernyak- 
Zhitnitsky DA we get 



<.^;'^"(Q^) = zV 



1 



2Q2 



91 



1 



2Q2 



ln(l + ^)-2 



(5.43) 



Despite its apparent complexity, this expression is very close numerically to the simplest interpolation 
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F;:'fJ{Q') 



1 



^/.(i + QV^' 



CZ\ 



(5.44) 



between the anomaly value at Q^ = and the pQCD result F^?^J^ '^ ^ {Q^) = | (47r/^/Q2) calculated for the CZ 
distribution amplitude. 

Thus, Eqs. (5.20), (5.44) model the modification of the basic /g-integral by power corrections. On the other 



hand, the modification of Iq by radiative corrections is described by Eqs.( 2.27| ),( ^.29 ). Though we obtained these 
two types of modifications in a completely independent way, it is tempting to combine them in a single expression. 
A self-consistent, but a rather time-consuming way to do this is to calculate the spectral density p'^""'''^(s, Q^) to 
two loops and apply the local duality prescription. Then both the radiative and power corrections would result 
from the same expression. We leave such a calculation for a future investigation. 

In the absence of a completely unified approach, we can try to get an interpolating formula by combining 
the two independent calculations described above. A natural idea is to write all the one-loop diagrams in 
the 6-representation a Id modified factorization and then substitute LpT^{x^\/h) by ^(a;,6) and the Born factor 
^Kii{\/SXbQ) by the modified version jS^Q^b^ K2{y/^bQ) . This will give a more reliable behavior in the small-Q^ 
region where the corrections are dominated by power terms. However, changing the structure of the Born factor 
would affect the radiative corrections and spoil the results at the high-Q^ end, where one should exactly reproduce 
the pQCD results. Since the perturbative corrections are rather small, we expect that a self-consistent inclusion of 
radiative corrections should be rather close to a simple product of the nonperturbative 1/(1 + Q'^ /a) factors and 
perturbative corrections from Eqs.(2.27),(2.29|). Such a product gives 



T^as 

7*7' 



AQ') 



1 



7rfn{l + Qyso) 



3 TT 



for the asymptotic form of the pion DA, and 

1 f 1 



f::4AQ') 



1- 



5a,(g2 



1 



1- 



lnQ2/A2 



(5.45) 



(5.46) 



1 



3n2 

5 



yso 



49 a,{Q^ 



108 TT j \\nQyAy 



lnQ2/A2^^"/8i" 



for the case when the pion DA (^^(a;;/i) coincides with ip'^ {x) for fi = Qo- These expressions have necessary 
interpolating properties: in the absence of radiative corrections they coincide with the local duality expressions, 
while for large Q^, when the power corrections can be ignored, they reproduce pQCD results. From Fig. 4, one 
can see that the curves for F"-^ j^oiQ^) and F^,^^a{Q'^) (with Qq « 0.5 GeV Q) in this model are sufRciently 
separated from each other which allows for an unambiguous experimental discrimination between them. 

It is instructive to make a more detailed comparison of the relative size of perturbative 0{as) and nonperturbative 
a/Q^ corrections. Taking A = 200 MeV, we observe that the perturbative correction for the asymptotic DA changes 
the lowest-order result by < 30% for Q^ > 0.5 GeV^^, This means that the pQCD expansion for the lowest-twist 
term in this case is self-consistent for Q^ as low as 0.5 GeV^. On the other hand, the power correction sq/Q^ 
exceeds 70% for all Q^ ^ IGeV^. This clearly indicates that pQCD results are not reliable below IGeV^. To 
reduce the ratio sq/Q"^ to the 20% level, one should take Q^ ^ 3GeV^. This is an illustration of the well-known 
statement (see, e.g., [p5[ ) that reliability of simplest pQCD formulas is limited in first place by power corrections 
rather than by the increasing value of the QCD running coupling as{Q^). The crucial fact here is that the scale 
So ~ 0.7 GeV^ determining the deviation from the pQCD 1/Q^ behavior is much larger than A^. It is also much 
larger than other typical nonperturbative scales like the square of the constituent quark mass M^ ~ 0.1 GeV^ or 

u2 \LD _ „„ /in ^ 0.07 GeV^). This observation 



the average transverse momentum {kj} (in the LD-model (5.25), {k 



so/10 ^ 



can be easily explained by the fact that k'j_ present in the modified Born term ( 5.28 ) is added to xQ^ rather than 
to Q^. This enhances the relative size of power corrections by a factor like l/{x). In full accordance with the 
statements made in refs. | |5^ , p3[ |, the onset of the Q^-region where the lowest-order pQCD result is reliable (in the 
sense that pQCD gives a good approximation) is determined by the size of the average virtuality xQ^ of the 
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FIG. 4. Combination \/2Q^_F^.^^o((3^)/47r (measured in GeV and equivalent to Q^F^,^^o{Q^), with the form factor 
F^,^^o{Q'^) normalized according to definition adopted in refs.[21,26,38]) as a function of Q^. The lower curve corresponds 
to our model with the asymptotic DA (Eq.(5.45)) and the upper one is based on Eq.(5.46). Data are taken from CELLO 
collaboration publication [14]. Preliminary CLEG data [15] (not shown) are very close to the lower curve. 

"hard" quark. If its value is too small, pQCD is unreliable even if the effective coupling as is negligible and 
perturbation theory for the lowest-twist contribution is self-consistent. 

VI. CONCLUSIONS 

In this paper, we discussed the status of QCD-based theoretical predictions for the F^*^^o{Q'^) form factor. As 
we repeatedly emphasized, in this case one deals with a rather favorable situation when QCD fixes both the Q^ = 
value (dictated by the axial anomaly) and the large-Q^ behavior governed by perturbative QCD. Still, constructing 
a dynamically supported interpolation between the two limits, it is very important to adequately reproduce at 
moderate Q^ the corrections to the asymptotic pQCD result, both perturbative and nonperturbative. 

Working within the framework of the standard pQCD factorization approach (SFA), which allows one to unam- 
bigously separate the contributions having different power-law behavior at large Q^ , we gave a detailed analysis of 
the one-loop coefficient function for the leading twist-2 contribution. To explore the role of the transverse degrees 
of freedom, we wrote the relevant Feynman integrals in the Sudakov representation and showed how the SFA 
produces the basic building blocks of the modified factorization approach (MFA) |23], such as the Sudakov- type 
double logarithms In (b) with respect to the impact parameter b± which is Fourier-conjugate to the transverse 
momentum k±. The fact that we derived the Sudakov effects within the lowest-twist contribution of the SFA, 
explicitly demonstrates that they should not be confused with the higher-twist effects. In other words, though 
the Sudakov terms are given by integrals over b± (or k±), they are purely perturbative and do not produce power 
corrections to the lowest-order pQCD result. 

Furthermore, we observed that the power corrections 1/Q^ due to the intrinsic transverse momentum are rather 
elusive both within the OPE-type factorization and the light-cone approach of Brodsky and Lepage. Contrary to 
naive parton expectations, the simplest handbag-type diagram in both cases does not produce an infinite tower of 
(l/Q^)" terms: such a series is generated by contributions corresponding to physical (transverse) gluons emitted 
from the hard propagator connecting the photon vertices. It goes without saying that an explicit summation of 
such terms is a formidable task in both of these approaches. A simpler picture emerges within the QCD sum rule 



30 



approach in which the infinite sum over the soft parts of the qG . . . Gq Fock components is dual to the qq states 
generated by the local axial current. An important observation establishing the connection between the QCD 
sum rule and light-cone approaches is that integrating the invariant mass s of the qq-pair over the pion duality 
interval < s < sq is. equivalent to using the effective two-body wave function ^^^(x, k±). The result obtained 
from the local quark-hadron duality (LD) ansatz applied to the lowest-order triangle diagram coincides with the 
Brodsky-Lepage interpolation formula 0], i.e., it reproduces both the Q^ — value specified by the axial anomaly 
and the high-Q^ pQCD behavior with the normalization corresponding to the asymptotic distribution amplitude 
for the pion. To test the sensitivity to the shape of the pion distribution amplitude, we proposed a model for the 
effective wave function '^^^{x, k±) which reduces to the desired DA after the fc_L-integration and still provides the 
correct limits for the form factor both at low and high Q^. 

In our analysis, the regions of small and large transverse momenta (responsible for power l/Q^ and Ug corrections, 
respectively) were studied separately, within the frameworks of two different approaches. In spite of this, the 
basic results written in terms of the fcj^-integrals look rather similar. A major challenge for a future study is the 
construction of a unified approach in which both the nonperturbative power-suppressed terms and the perturbative 
radiative corrections emerge from the expansion of the same expression. The quark-hadron duality approach 
provides a framework in which such a self-consistent unification is guaranteed. The only missing ingredient is the 
perturbative spectral density p'^'"°'^'^{s,Q^) at the two-loop level. 

There are two further improvements which should be made in the perturbative part of the problem. First, it is 
necessary to fix the argument of the running coupling constant ag. In our analysis, we cither left it unspecified and 
estimated the corrections assuming that as/ir « 0.1 or took A = 200 MeV in the 1-loop expression for as{Q'^). 
However, for a precise comparison with experimental data, estimating the magnitude of the ag-correction one 
should explicitly specify the UV renormalization scheme, fix the parameter ^n in the argument of the running 
coupling asinu) and use the proper value of the QCD scale A. A very effective scale-fixing prescription is provided 
by the Brodsky-Lepage-Mackenzie approach |7Q]. To use the BLM prescription, one should calculate two- loop 
pQCD corrections to the coefficient function containing quark loop insertions into the gluon propagator. Another 
problem is the inclusion of the effects due to the two-loop evolution of the pion distribution amplitude 71-73|. 



Originally, the relevant corrections expanded in terms of a few lowest eigenfunctions of the one-loop kernel, were 
found to be tiny [g[ . A recent progress [Q in understanding the structure of the two- loop evolution suggests that 
higher harmonics cannot be neglected, and the size of the two-loop evolution corrections is somewhat larger than 



estimated in B. However, our preliminary numerical estimates 75 1 of the effects due to the modified evolution 



developed in ref. [[76| do not indicate appreciable changes for the /-integral. 
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